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Abstract 

Few layer graphene systems such as Bernal stacked bilayer and rhombohedral (ABC-) stacked 
trilayer offer the unique possibility to open an electric field tunable energy gap. To date, this energy 
gap has been experimentally confirmed in optical spectroscopy. Here we report the first direct 
observation of the electric field tunable energy gap in electronic transport experiments on doubly 
gated suspended ABC-trilayer graphene. From a systematic study of the non-linearities in current 
versus voltage characteristics and the temperature dependence of the conductivity we demonstrate 
that thermally activated transport over the energy-gap dominates the electrical response of these 
transistors. The estimated values for energy gap from the temperature dependence and from 
the current voltage characteristics follow the theoretically expected electric held dependence with 
critical exponent 3/2. These experiments indicate that high quality few-layer graphene are suitable 
candidates for exploring novel tunable THz light sources and detectors. 
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Few layer graphene (FLG) -stacked sheets of carbon atoms on a honeycomb lattice- are 
materials with unprecedented properties such as an electric held tunable band structure 
which depends on the stacking order of the layers Bernal (AB-) stacked bilayers and 

rhombohedral (ABC-) stacked trilayers open a band-gap when the energetical symmetry 
between the outer layers of the FLGs is broken [MB [I6H2H], e.g. by an external electric 
held (see Fig. & and b). Furthermore, ABC-stacked trilayers are also an ideal system to 
study competing physical aspects of the single-particle and many-body physics, with the 
observation of rich phase diagrams [IBII3]. 

The unique possibility to continuously tune the band structure in FLGs simply by means 
of gate voltages is a novel functionality which is attracting a growing interest in fundamental 
science and in future graphene-based applications [22]. In particular, a gate-tunable band- 
gap, characteristic of few layer graphene materials, is at the core of optoelectronics. So far, 
only infrared spectroscopy experiments have directly measured the value of this electric held 
induced band-gap in both AB-bilayers [I9H2I] and ABG-trilayers [10]. These experimental 
hndings are contrasted by the lack of the direct observation of a tunable band-gap in charge 
transport measurements where a transport gap at a very diherent energy scale has been 
reported [231 - 126] . in spite of the observation of an electric held induced insulating state 
indicating the opening of an energy gap [2Z]- Proper consideration to the role of disorder 
induced sub-gap states can shed light onto this dichotomy between optical spectroscopy and 
electrical transport reports. In particular, the infrared aborption is largely dominated by the 
band-to-band transitions rather than transitions from the small density of disorder induced 
states. On the other hand, in electrical transport the hopping of electrons between sub¬ 
gap impurity states dominates the current how hindering the observation of the intrinsic 
energy gap. Indeed, the characteristic energy scale (activation energy) associated to the 
hopping conduction mechanism is related to the ionization energy of the sub-gap states 
[29] and brings no knowledge of the intrinsic band-gap opened in the material. The direct 
observation of the FLG gate tunable band-gap in electrical transport experiments is still an 
open quest which is fuelling theoretical discussions [30] . 

Here we report the direct observation of an electric held induced tunable band-gap in 
ABG-stacked trilayer graphene with electrical transport measurements. A systematic study 
of the temperature dependence of the minimum conductivity in these devices reveals that 
conduction is due to thermally excited charge carriers over the energy gap and not to hopping 
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between sub-gap impurity states. We find that the magnitude of the measured band-gap 
depends on the intensity of the external perpendicular electric displacement with values of 
4.6 meV for an average electric displacement of -110 mV/nm. The values of the energy gap 
estimated from the temperature dependence of the minimum conductivity correlate well 
to the energy range of non-linearity measured in the current-voltage (I-V) characteristics. 
These observations stem from the high quality of our devices which results in the absence 
of any hopping mechanism assisting the electrical conduction. 

ABC-stacked trilayer graphene flakes are obtained by micromechanical exfoliation of pris¬ 
tine natural graphite onto Si02 (295 nm)/p-Si. The number of layers and the stacking order 
are identified by Raman spectroscopy [mEHEg. In particular, the 2D-peak of ABC-stacked 
trilayers is asymmetric as shown by the relative heights of the 6 Lorentzians fitted to it (see 
Fig. IJ;). Source and drain contacts are patterned by standard electron-beam lithography, 
metal deposition of Cr/Au (lOnm/lOOnm) and lift-off process. The independent control 
of the external perpendicular electric field acting onto the trilayer, which determines the 
opening of a band-gap (see Fig. R). and of the charge density is readily obtained in double 
gated devices by changing independently the two gate voltages |22] . To achieve high quality 
devices, suspended double gated structures are fabricated [33] (see Fig. W-f and supple¬ 
mentary information). The width and length of the top-gate electrode are much larger 
than the graphene channel length and width, therefore ensuring that a uniform electric 
field acts on the conductive channel. Typically, the charge carrier mobility in these devices 
is /i > 25000cm^/Rs at 4.2K and at a charge concentration of ~ 5 * with no 

noticeable residual doping. 

Electrical measurements were performed in a He^ refrigerator in the temperature range 
between 300 mK and 100 K. The conductance of the suspended double gated devices was 
measured using a lock-in amplifier in a voltage-biased two-terminal configuration. The choice 
of the excitation voltage is pivotal to understand if an insulating state with an energy gap 
smaller than the recently reported values due to many-body states [HEI is present in ABC- 
trilayer. Hence we have progressively lowered the values of excitation voltage down to the 
lowest possible values which did not affect the reading of the maximum of resistance at any 
of the measured temperatures. In this way, we have also prevented heating of the electrons 
and the occurrence of non-equilibrium effects which could easily hinder the observation of 
small values of the intrinsic energy gap. 
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FIG. 1. (a) Schematic crystal structure of rhombohedral (ABC-) stacked trilayer graphene (left) 

with the first neighbour in-plane 70 and out of plane 71 tight binding parameters highlighted, (b) 
Energy dispersion for ABC-trilayer with zero (left) and finite (right) external perpendicular electric 
field, (c) 2D-Raman peak for ABC-trilayer measured with a 532nmn laser, 5mW power and a spot 
size of 1.5 fim. The dots are the experimental data points, whereas the red continuous line is a fit 
to 6 Lorentzians (continuous blue lines), (d) Schematic cross section of a suspended double gated 
device and the electrical measurement configuration. Panels (e) and (f) show false color scanning 
electron micrographs of suspended double gated graphene devices, and a test suspended back-gated 
graphene structure fabricated on the same chip. 
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FIG. 2. (a) Zero-bias resistance vs. Vhg measured at T = 0.3 K and for different fixed Vtg as 

indicated in the graph, (b) 3D Plot of R vs. Vbg and Vtg. The dashed line is a guide line for the 
eyes connecting the maximum of resistance Rmax- (c) R vs. Vbg for Vtg = 15V (left) and Vtg = OF 
(right) and for different temperatures. 


Figure]^ and b show the resistance of a snspended donble gated ABC-trilayer (channel 
width and length of 500nm) measnred at 300 mK as a fnnction of back-gate voltage (Vbg) 
for different fixed top-gate voltages (Vtg). We define the average external perpendicnlar 
electric displacement as D = (Dbg + Dtg)/2 with Dbgda + Dbgdsi 02 /^ = Hg and Dtgdtg = Vtg 
(dtg = ISOnm, da = S5nm, dsi 02 = 215nm and £ = 3.9 for Si02). Upon increasing D, 
the maximnm of resistance (Rmax) increases from 7.9 KOhm np to 1 MOhm and shifts its 
position in the gate-voltages, see Fig. If. This insnlating state is at least ten times more 
resistive than the recently reported insnlating state dne to many-body effects [niE] and its 
dependence on the voltage applied to the two gates is snmmarized in the color plot of Fig. 

We find that Rmax moves linearly as a fnnction of Vbg and Vtg, as expected for constant 
capacitive conpling to the two gates. At the same time, the temperatnre dependence of R 


5 























FIG. 3. (a) Color plot of log(i2) as a function of Vtg and Vhg. The inset shows I-V characteristics at 

the Rmax for 3 different gate configurations highlighted in the main graph. I, II and III correspond 
to 0.05, 0.08 and 0.13 V/nm electric displacement. The dashed lines in the inset highlight the 
non-linear energy range assigned to the intrinsic gap (A=3mV for D=0.08 V/nm). The differential 
conductance dl/dV versus Vsd and Vhg for fixed Vtg corresponding to I, II and III is plotted in 
panel (b). The sketches on the right hand side show a schematic low-energy band structure of 
ABC-trilayer and the position of the Fermi level [Ep) for the three different values of I4g indicated 
by the white dashed lines. 


versus Vtg shows that the larger is D the more pronounced is the temperature dependence 
of Rmax, see Fig. HI- Both the increase of Rmax as a function of D and its pronounced 
temperature dependence at large D suggest the opening of an electric-held induced band- 
gap due to the breaking of the inversion symmetry between the outer layers of ABC-trilayer 
graphene [iHin]. 

Another indication of an electric held induced band-gap comes from the current voltage 
characteristics. The inset in Fig. shows I-V curves for three different perpendicular 
electric helds conhgurations highlighted in the color plot. In all cases, the I-V curves display 
an insulating state at low Vgd bias, though for very different AVgd ranges. In particular, AVgd 
decreases with decreasing D, i.e. AVj^ < AVJJ < AVJJ^. Correspondingly, the differential 
conductance (dl/dV) measured as a function of Vgd for the same gate conhgurations is low 
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when the Fermi level is in the middle of the band-gap. In contrast, a high dl/dV is observed 


when the Fermi level is in the conduction or valence bands, see Fig. 

To demonstrate that we probe an intrinsic band-gap in these electrical transport exper¬ 
iments, we performed a detailed study of the temperature dependence of the conductivity. 
These measurements can distinguish if charge carriers are hopping between sub-gap impu¬ 
rity states or if they are thermally excited between valence and conduction band over the 
energy gap (A). In the hopping regime, depending on the specihc conduction mechanism 
and on the relevance of correlation effects, the temperature dependence of conductivity (a) 
for a two dimensional system is described by a = cToexp(—(T q/T)") with n = 1, 1/3 or 1/2 
corresponding to nearest neighbour hopping, two-dimensional Mott variable range hopping 
and Efros-Shklovskii variable range hopping in the presence of Coulomb interaction between 
the localised states. Tq is a characteristic parameter which is a function of the activation 
energy of hopping izg. On the other hand, in the absence of impurity states thermally 
activated electrical transport described by cr = cToexp(—A/2fcsT) (fcg Boltzman constant) 
is expected. In addition, for short channel high quality devices and at low temperatures 
{ksT < A), a sub-gap electrical current is expected due to quantum tunneling through the 
energy barrier dehned by A. 

Figure]^ shows some representative Arrhenius plots of the temperature dependence of the 
minimum conductivity a for different values of D. It is apparent that in all cases the exper¬ 
imental data are well described by a linear £t, as expected for thermally activated transport 
in intrinsic semiconductors and A is the only htting parameter for the slope. We hnd that 
the band-gap can be continuously tuned up to A = 4.6±0.06mel/ at D = —llOmV/nm, see 
Fig. 4b. To understand the temperature dependence of a down to cryogenic temperatures 
(see Fig. 4c), we consider the combination of thermally activated transport and quantum 
tunneling through a short-channel semiconductor with a clean band-gap. The conductivity 
is given by the relation: 



( 1 ) 


A/2 


where DOS{E) is the density of states for ABC-trilayer, f{E) is the Fermi distribution 
and cTi is the measured low-temperature residual conductivity due to quantum tunneling (see 
supplementary information). The low-energy dispersion of ABC-trilayer is described by one 
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FIG. 4. (a) Arrhenius plot for different values of external electric displacement. The dots are 

experimental data and the continuous lines are fit to thermally activated transport over the band- 
gap A. e on the conductivity axis is the constant of Napier (b) Shows the estimated values of 
the band-gap as a function of D. The continuous line is a ht to A oc (c) Plot of the 

minimum of conductivity as a function of inverse temperature for different average external electric 
displacements, as indicated. The dots are experimental data points, whereas the continuous line is 
a ht to Eq. 1. (d) Linear ht of /n(A) vs ln{D) with best ht A = 3 ± 0.2. e on the axis of A is the 
constant of Napier. 

pair of conduction and valence bands trigonally warped for E < lOmeV evolving into a cubic 
dispersion [E oc k^) for E > lOmeV^ [IHTj. At E = lOmeld the Fermi surface undergoes 
a topological discontinuity known as Lifshitz transition HEIEIDS]. In our experiments, 
both the temperature range and the values of A (see Fig. 4b) are smaller than the energy 
corresponding to the Lifshitz transition {El). Therefore, we are probing electrical transport 
at the edges of the conduction and valence bands. The opening of a band-gap introduces 
a parabolicity at the band edges for each separate pocket known to exist in the low energy 


















dispersion of ABC-trilayer |lj, resulting in a constant density of state equal to m*/2nh‘^ with 
m* effective mass for each pocket. Figure]^ shows that considering a constant DOS{E) in 
Eq. (1) gives an accurate £t to the experimental data using as the only fitting parameter 
A since the mobility is obtained from the held effect measurements, the effective mass 
{m* = 0.05 ± O.Olmo, with mo free electron mass) is kept hxed after the hrst £t to the data 
at D=100 mV/nm and cxi is the low temperature residual conductivity. For A El the 
low energy dispersion is similar to 3 gapped Dirac cones (i.e. DOS{E) \E\ ■ 6{\E\ — A|)). 
However, even in this limit, the difference between the values of A found considering a linear 
or parabolic dispersion is within experimental error. The estimated values of A from the £t 
to Eq. (1) are consistent with the values found from the Arrhenius plot analysis and they 
also correspond to the energy scale of the non-linearity measured in the I-V characteristics 
(see Fig. 1^ and Fig. |^). The self-consistent tight binding model predicts that the low- 
electric field dependence of the energy gap in ABC-multilayers is A oc with N the 

number of layers 0171 . We verify this theoretical prediction by htting the experimental data 
for ln{A) vs. ln{D) to a linear slope. The best fit gives N = 3 ± 0.2 as expected for the 
intrinsic band-gap in ABC-trilayer graphene, see Fig. B- 

Finally, we note that both Nearest Neighbour Hopping (NNH) and thermally activated 
transport have a similar exponential decay of the resistivity as a function of temperature 
[2ni El] • However, in our devices the electron mean free path at the edge of the conduction 
and/or valence bands is typically > 200nm -comparable to the transistor channel length 
and width (both ~ 500nm). Therefore, the electrical transport regime in our experiments 
is quasiballistic and not diffusive as in NNH. If we assume that A is the activation energy 
of NNH, we hnd a localization radius of lOnm, i.e. one order of magnitude smaller than 
the electron mean free path in these devices. Consistently, we find that the values of A 
estimated from cr(T) correspond to the energy scale of the non-linearity measured in the 
I-V characteristics, see inset in Fig. 1^ and Fig. i>. contrary to NNH regime PI EU. 
Furthermore, we estimate the same values of A in devices with different aspect ratios (see 
supplementary information), conclusively demonstrating that our experiments provide a 
direct measurement of the gate-tunable band-gap in transport experiments in ABC-trilayer. 

In conclusion, we have reported a detailed study of the source-drain bias, temperature and 
perpendicular electric held dependence of the resistance in suspended double gated ABC- 
trilayer graphene. All these independent measurements consistently show the appearance of 
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an energy gap in the density of states of ABC-trilayer. We demonstrate that we can tnne 
A np to 4.6mef^ for an average electrical displacement of llOmld/nm. Onr experiments 
present the hrst direct observation for a clean gate-tnnable band-gap in electrical transport 
measurements in ABC-trilayer graphene. These hndings indicate that high quality few-layer 
graphene are an attractive system for the developement of gate tunable and highly efficient 
THz sources and detectors where a clean energy gap is an essential requirement. 
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Supporting information. Additional information on fabrication of suspended and 
doubly gated ABC-stacked trilayer graphene, current annealing of suspended and doubly 
gated ABC-stacked trilayer graphene, energy gap in other doubly gated ABC-stacked tri¬ 
layer graphene devices and low temperature quantum tunnel conductivity is available in 
the Supporting Information. This material is available free of charge via the Internet at 
http://pubs.acs.org,. 
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